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SUMMARY 

A theoretical and experimental study of slender tubular columns for 
possible use In space structures Is conducted In the presence of partial rotah- 
tlonal end restraints. New explicit formulas are derived for computing the 
buckling load and the lowest natural frequency of perfectly straight uniform 
elastic members with rotational end restraints possessing linear moment- 
rotation characteristics. An exact solution In the form of a transcendental 
equation, and a numerical solution using second-order finite-differences are 
also presented. The presence of an Initial Imperfection Is also Incorporated 
Into the numerical procedure. Each of the solutions presented originates from 
the partial differential equation and the associated boundary conditions for 
the problem. It Is shown that the buckling load and the natural frequency 
formulas developed are highly accurate. Vibration tests are conducted on an 
Icperfect tubular steel member In the absence of an axial load. A darling 
concept consisting of a string-mass assembly Is explored In a few of the tests. 
In addition to the structural danping present In the system. Three passive 
damping configurations were Investigated; one with three small equidistant lead 
shots attached to a string Inside the tubular member, one with the three lead 
shots attached at the center of the tube, and the other with one lead shot at 
the center of the tube. The three lead shot configurations both provided con- 
siderably greater damping than the single lead shot damper. 
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SYMBOLS 

damping coefficient 

outer diameter of the tubular member 

Young's modulus 

lowest natural frequency of the column 
experimental natural frequency 
dimensionless natural frequency 
panel length 
moment of inertia 
end rotational stiffness 
dimensionless end rotational stiffness 
n»rober length 
axial load 

buckling load of the column 
dimensionless buckling load of the column 
time 

wall thickness of the tubular menber 

net lateral displacement 

initial out-of-straightness 

member longitudinal ordinate 

axis parallel to the face of the proximity probe 

axis perpendicular to the face of the proximity probe 

envelope of the raidspan maximum dynamic deflection 

time interval 

danping ratio 

mass per unit length of the column 
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The analysis for buckling and natural vibration of elastic coluons with 
pinned or fixed end conditions Is well-known (see refs. 1 and 2). The prac- 
tical columns In real space structures, however, seldom possess such Ideal 
boundary conditions. The connections and the adjoining members which frame In 
at the ends of a typical column In a space structure provide rotational 
restraints which affect the behavior of the column Itself. The buckling load 
of a partially restrained column may be obtained by using the classical eigen- 
value approach leading to a transcendental equation (see for exaaple, ref. 3) 
which mist be solved Iteratively. Also, It Is possible to formulate an exact 
solution for determining the natural frequency of such a member In the form of 
a complicated transcendental equation which must also be solved Iteratively. 
This paper presents an outcome of a theoretical and experimental Investigation 
of tubular members restrained partially. New explicit formulas are derived for 
conputing the buckling load and the natural frequency of perfectly straight 
uniform elastic members having partial rotational end restraints possessing 
linear moment-rotation characteristics. An exact solution In the form of a 
transcendental equation, and a numerical solution using second-order finite- 
differences are also presented. The presence of initial imperfection is also 
incorporated into the numerical procedure. Each of the solutions presented 
originates from the partial differential equation of motion governing the 
dynamic response of the column in the presence of a constant axial load and 
the boundary conditions representing the rotational springs. The vibration 
tests are conducted on an inperfect tubular steel member with rotational 
elastic end restraints and in the absence of an axial load. The use of a 
damping concept is explored in a few of the tests, in addition to the 
structural danping already present in the system. 
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PROBLEM STATEMENT AND GOVERNING EQUATIONS 
Figure 1 shows schematically a slender column of length L with a hollow 
circular cross section. The outer diameter of tlie section Is Dq and the wall 
thickness Is tg. The column has an Initial out-of-straightness given by w. 

A static axial load P Is applied to the column. Also, the ends of the co1i«m 
are not allowed to translate laterally, and possess partial rotational 
restraints each having a rotational spring stiffness of K ir»-lb/ radian. The 
material of the coliam Is elastic. The governing differential equation of 
equillbrijn for the member shown In figure 1 may be written as: 


ax 
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a w 

2 

ax 


2 - 

a w 


ax 


at' 


+ c 


at 


( 1 ) 


In which 

w(w,t) = lateral displacement 

w(x) = Initial out-of-straightness, or Imperfection 

El = flexural rigidity of the column 

P = axial load 

p = mass per unit length of the column 

C = damping coefficient 

Equation (1) without the terra Involving w can be found In the standard 
texts on structural dynanlcs (see, for example, ref. 2). As evident from the 
last term In equation (1), the system dancing Is assumed to be proportional to 
the velocity at ar\y given time. 

BOUNDARY AND INITIAL CONDITIONS 

The boundary conditions for the problem shown In figure 1 are as follows: 


w(0,t) » 0 
w(L,t) “ 0 


( 2 ) 

( 3 ) 
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El ~ (O.t) » K ^ (O.t) ( 4 ) 

ax 

2 

El (L,t) “ -K — (L,t) (5) 

Equations (2) and (3) represent zero lateral translation of the aeaber ends. 
Equations (4) and (5) are the natural boundary conditions dependent upon the 
flexural stiffness of the column ends and the restraint aocaents developed by 
the rotational springs. 

The initial conditions for obtaining the natural frequency of the column 
are taken as follotrs: 


w(x,0l = <^(x; P. K. El, 1 ) (6) 

U (x.O) = 0 (7) 

in which ^ is a function of x and the terms P, K, El, and L. The 
initial condition given by equation (6) states that at tine t equal to zero, 
the column has been displaced ("plucked”) resulting in a shape given by a func- 
tion of X, the applied axial load P, the spring constant K, the flexural 
rigidity El, and the length L. Equation (7) imposes a zero initial velocity 
condition. 


EXACT SOLUTION 

For the case of no initial imperfection, that is, w “ 0. equation (1) can 
be solved by the separation of variables technique by letting; 

w(x,t) " W(x) T(t) (8) 

Substituting equation (8) into (1) leads to the following pair of ordinary 
differential equations (see ref. 4): 
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d^W 2 d^W 

■^ + k - XW - 0 

^ dx^ 


d^T C dT 2 ^ 

dt2 p dt 


(?) 


( 10 ) 


In which X Is a parameter to be determined, u Is the undamped circular 
frequency given by: 


iji 



EIX 


and k Is given by: 

I El 

The exact solution to equation ( 9 ) Is given by (ref. 4 ): 


( 11 ) 


( 12 ) 


W(w) = cos + B2 sin pjX + Bj cosh 5lnh 


1n which 


2 ^ k^ ±/k^ + 

P1.2 ■■ * r- 


4 x 


( 13 ) 


(14) 


and B^ through B4 are constants of Integration. The solutions to equa- 
tion (10) are all of the form (ref. 6): 


T = F(t) 


(15) 


In which F(t) depends upon the problem Input parameters, u, and the Initial 
conditions (6) and ( 7 ) Invoked via equation (8). The damped circular frequency, 
uj, Is given by (ref. 6): 


Uj ° 0^1 - c 


( 16 ) 


1n which c Is the damping ratio. 
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The solutions based on equation (13) for Idealized pinned or fixed end 
conditions are given In the literature (see ref. 4). Here, an exact solution 
for a column with partial restraints and having the boundary conditions given 
by equations (2) through (5) Is presented. Using equation (8), these condi- 
tions reduce to the following ones: 


W(0) = 0 


(17) 


W(L) “ 0 


(18) 


d^W 


El ^ (0) n K (0) 


(19) 


d^W 


El ^ (L) » -K (L) 


( 20 ) 


Substituting equation (13) Into equations (17) through (20) results In a matrix 
eqjatlon of the type: 


^11 

ai2 

«13 

*14 




^ N 

0 

821 

«22 

^23 

*24 

< 

82 


° 

®31 

«32 

*33 

*34 


S3 

0 

_«41 

842 

*43 

*44_ 




.0. 


In which the a^j terms are given In appendix A. Setting the determinant of 
the matrix In equation (21) equal to zero leads to a complicated transcendental 
equation. For a given value of the axial load P, the colixnn dimensions and 
material properties, and the spring constant K, the resulting transcendental 
equation may be solved Iteratively for the lowest nontrivial value of x. A 
substitution of this x value Into equation (11) gives the circular frequency 
of the partially restrained column. To find the coluun static buckling load, 
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however, the smallest nontrivial value of P oust be sought with X ■ 0, since 
the frequency of a colictn becomes zero when the axial static load equals the 
buckling load. 

When the axial load P Is less than the buckling load, the deflection 
versus time relationship may be found by substituting the u> value, ascer- 
tained by the iterative procedure, Into equations (13) and (15) which. In turn, 
must be substituted Into equation (8). 

The disadvantage of the exact solution Is that It can be obtained only 
through trial and error. In the following section, explicit fortaulas for 
computing the buckling load and the natural frequency of the partially 
restrained column are derived. 

BUCKLING LOAD AND NATURAL FREQUENCY FORMULAS 

In this section, new explicit formulas for calculating the buckling load 
and the natural frequency of the partially restrained colmn shown In figure 1, 
In the absence of Initial liqjerfectlons, are derived and are based on a highly 
accurate assuned mode shape of the column. The procedure Involves selecting a 
trigonometric expression which satisfies the boundary conditions a priori, and 
then using Galerkins' orthogonality criterion. 

Let the solution to equation (9) be taken In the following fora: 

W = s1np+ A2 (i - cos-?p) (22) 

In which A^ and A 2 are constants. This expression satisfies the boundary 
conditions (17) and (18). The sine term alone would represent a pinned meaber 
while the term In the parentheses approximates a completely fixed end oeaber. 

A combination of the two functions has not been explored by any previous Inves- 
tigator for simulating the deflected shape of the member with end conditions 
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somewhere between pinned and fixed cases, that Is, for the partially restrained 
end conditions. Equation (22) represents the approximate fundamental node of 
the member, thus the buckling load and the natural frequency expressions are 
associated with this mode. Furthermore, due to the symmetry of the W 
expression about the member midspan. It Is sufficient to enforce only one of 
the conditions given by equations (19) and (20). Substituting W from 
equation (22) Into equation (19) results In the following relationship 
between Aj^ and A 2 : 


A A 

A o ^ — A 

2 4irEI 1 


Substituting equation (23) Into (22) gives; 


(23) 


W = 


. r j ’tX KL f . 2tix\ 


( 24 ) 


If equation (24) were exact solution, It would satisfy equation (9) exactly. 
Since It 1s an approximate solution, Its substitution Into equation (9) results 
In an error expression given by e as follows: 


\|^^(s1n — - 16n cos _j + — ^-sln— + 4n cos-j-j 
sin ^ + n^l - cos ^ 


(25) 


In which: 


KL 

“ 4^EI 

Invoking use of Galerkin's orthogonality criterion gives the following 
equation: 


(26) 


/ 1“ 


dx o 0 


(27) 


0 


5- K-. 
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Upon performing the Integration, rearranging, and noting that Aj y 0 for the 
nontrivial solution, an explicit expression for X Is obtained from equa- 
tion (27). Using this X expression, and noting that u “ 2tf, where f Is 
the undamped natural frequency In Hertz, the following formula Is obtained: 


f 


1 

El 

- S* '' '"''(I ■ r)* • n) 


P 

12(»EI)^ + 32 EIKL + 7 OCL)^ 

•* 


(28) 


Of course, the damped frequency can be obtained by using equation (16). 

A buckling load formula for the column shown In figure 1 can be obtained 
by setting f from equation (28) equal to zero, since the frequency of the 
column approaches zero as the axial load ap, *oaches the buckling load. Thus, 
setting the right side of equation (28) to zero, simplifying, and solving for 
the critical P value leads to; 


P 

cr 



12(ttEI)^ 80 EIKL 12(KL)^ 

12(TrEI)^ + 32 EIKL + 3(KL)^ . 


(29) 


Formulas (28) and (29) can be used to predict the natural frequency, and the 
elastic buckling load of a perfect column with equal partial end restraints. 

The formulas derived above are directly applicable In the absence of an 
Initial lo^erfectlon. It Is well-known (ref. 3) that the maximum load which an 
Imperfect elastic column may carry. In the oreaence of pinned or fixed end 
conditions, asymptotically approaches the Euler buckling load of the corre- 
sponding perfect roen*er. Although the details are not given here, It can be 
shown that the same conclusion Is applicable to the partially restrained coliaan 
considered herein. Furthermore, similar reasoning Is applicable to the evalua- 
tion of the lowest natural frequency of both perfect and Imperfect members. 
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COW>ARISON FOR LIMITING BOUND-IRY CONDITIONS 


If the end spring stiffness K Is taken as zero, the limiting case of a 
column with pinned ends Is obtained. It can be verified that If K Is set 
equal to zero In equations (28) and (29), the resulting f and values 
are the exact solution for pinned ends (see refs. 1, 2, and 5). For a fixed 
end column, the value of K equals Infinity. To evaluate f for this Halting 
case, the numerator and the denomi nator terras under the radical sign should 
first be divided by and the resulting expression evaluated as K 
approaches Infinity. Using the expression for f thus obtained, the circular 
frequency u for a member with zero axial load Is found to be as follows: 


u = 22.79 


\l 


El 

pL 


(30) 


which Is only slightly different from the exact u expression given In 
reference 5 with a coefficient of 22.40 outside of the radical te To 
evaluate P^p with an Infinite value for K, the numerator and the denocslnator 
terms on the right side of equation (29) should first be divided by and 
the resulting expression evaluated as K approaches Infinity. This process 
leads to the exact buckling load of a fixed end column given In reference 1. 


INTERACTION RELATIONS 

Figure 2 presents nondlraenslonal Interaction curves relating the nondlraen- 

- 2 

slonal frequency f and Its square (f] with no axial load (P » 0 In 
fig. 1), and the buckling load P , to the nondlraenslonal spring constant K , 
defined as follows: 

(31) 



00 


( 32 ) 
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in which f Is given by equation (28) with P “ 0; f„ Is the value of f 

with K ° Is given by equation (29); and P„ Is the value of P^^ 

with K “ As seen from this figure, f Increases from about 0.433 tc 1.0 

- 2 

as K Is varied between zero and Infinity. The variation of (f) with K 

Is also shown. Similarly, P Increases from 0.25 to 1.0 as K Is varied 

between zero and Infinity, viz., from the pinned to the completely fixed end 

- 2 

condition. Substantial Increases In f, [f] , and P occur as K Is 

increased from zero to about 20, beyond which only small Increases are observed 

even for very significant Increments In K . This shows that what mAy appear 

to be a nominal end restraint may actually provide a considerable degree of end 

fixity rather than a nearly pinned condition. In fact, a 'moderate' amount of 

end restraint may provide near-complete end fixity. 

As may be seen from equation (28), the presence of an axial load P In 

the presence of partial end restraints would reduce the natural frequency f. 

- 2 

Thus, the f, and (f) versus K relations with P > 0 would fall below 
those shown In figure 2. 

A comparison of the P versus K relation shown In figure 2 was made 
with a finite element solution of the problem where the column was divided Into 
a total of ten segments. The results were found to be almost Identical. The 
validity of the formula (29) was also verified by confuting a few of the points 
on the P versus K relation In figure 2 using an exact eigenvalue solution 
discussed earlier In this paper. Again, the results were almost Identical. 

FINITE-DIFFERENCE SOLUTION 

Using second-order central finite-difference expressions (see ref. 7), the 
partial differential equation (1) may be written as follows; 
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El 


4 ^ V2,j^ 


.2 J ^ VlJ^ * 

n 


(w, ^ - 2w, , + w, ♦ 


(At) 


2 ‘"ij-l ""ij "i,j*l 


(At) 


— (-W 4' W ) ° 0 

2 1J i.J+r 


(34) 


In which: 

h ° panel length along the x-axIs of the coltan (fig. 1) 

At = tlw Interval 

4.(x^) ° second derivative relative to x of tl« 1nq>erfect1on 
function w 

x^ ° 1h, for each 1 = 1, 2, 3, ... 

The subscript 1 refers to the 1th panel point over the dowln 0 < x < L, and 
the subscript j refers to the miaber of tiae Increaents such that the tiae 
at J Is given by the following equation: 


tj ° j(At), for each j = 0, 1, 2, 3, ... 

Slollarly, the boundary conditions (2) through (5) can be expressed In the 
f 1 n1 te-dl f ference fora as follows: 


( 

( 


w_ . = 0 
O.j 

(35) 

w 0 

MJ 

(36) 

E* ^ ^ /El K\ 

— + vlW. .+(r: TTiW, .*’0 

(37) 

h z) -IJ Vh Z) l,j 


El K\ /El K\ 


h 2 ) Vl J 2 ) VlJ “ ° 

(38) 


Applying equation (34) at 1 “ 1, 2, 3, .... (M-1) and Invoking the condi- 
tions (35) through (38) leads to the following oatrlx equation: 
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1 

c ® - 

1 b, + b. 

(40) 

3 4 


C2 » bjci 

(41) 

^3 “7^ 

(42) 

(At) 


b. 

4 At 

(43) 

p^ “ PCj4,(x^) 

(44) 


and [r] is a syaaetric coefficient matrix of the order (M-1) by (M-1). 
Appendix B defines the terms of this matrix for H => 6 as an example. 

Equation (39) m^y be used to predict the lateral deflections 
j and w-i j_i are know. To avoid having a negative time interval due to 
the use of the central finite-difference scheme i^hen starting the recursion 
indicated in this equation, the following special forward start-up differertce 
equation for is derived in appendix C: 



•“eVz * ('^6 * * (^"6 ■ 

* (-«>6 * 


in which: 


b 


6 


o 



El 

P 


b 


7 



P 


(At)^ 

n . ----- - 

8 2 


(45) 


(46) 


(47) 


( 48 ) 
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In using the recursive equation (39), the panel length h and the tine incre- 
ment At should be kept small to avoid numerical instability of the solution. 
For the member tested in the laboratory described later in this paper, it was 
found sufficient to use h = L/6, and At = 0.001 sec to ensure both con- 
vergence and numerical stability. 

EXPERIMENTAL INVESTIGATION 

The experimental portion of the investigation consisted of conducting 
several natural vibration tests on a tubular steel member of about 12 ft 
length, and outer and inner diameters of 0.50 and 0.37 in., respectively. The 
ends were partially restrained by specially designed connections, and no axial 
load was applied. The steel mend>er was used since very slender graphite com- 
posite tubes are not yet available, and the procedures used in this investiga- 
tion are good for both materials. The ffrtanber had some initial imperfection due 
to manufacturing and handling processes. 

Figure 3 shows sketches of the mender tested, and the details of the end 
connections. Each end connection was a two-piece assembly. One piece was a 
steel blade 1-1/2 in. long with a 3/16 in. x 1/2 in cross section welded to a 
steel base plate. The other piece was a steel clevis designed to fit snugly 
over all but 1/4 in. of the blade and fastened to it by screws fitted through 
two drilled holes. The exposed 1/4 in. of the blade, thin in one direction 
(3/16 in.) but thick in the other (1/2 in.) provided partial rotational 
restraint in one plane of vibration and almost total rotational restraint in 
the orthogonal plane. This design provided the member a preferred direction of 
vibration, namely, about the weaker axis of the blade. The top end of the 
clevis was a 1-in. threaded rod which screwed into the steel raejuber. The 
member and the end fixtures were threaded so that one end had right-hand 
threads and the other had left-hand threads. These opposite thread directions 
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allowed the systea to act like a tumbuckle so that the top and the bottoa end 
fixtures could be easily aligned during assembly. The upper base plate was 
bolted to a heavy bracket which In turn ms bolted to a backstop, and the lower 
base plate was bolted to the floor directly beneath the upper one. The clevis 
portions of both end fixtures were screwed Into the appropriate ends of the 
member. The upper end of the member was then attached by the clevis-blade 
Joint to the base plate. The turnbuckie feature of the system was then used to 
align the attachment holes on the lower end clevis-blade Joint so that no axial 
tension or coi^resslon was Induced In the systea. 

The Instnoentatlon used In the tests consisted of a proxlolty probe, 
proximity probe target, and a strip chart. The proxlnl^ probe was mounted on 
a bracket at midspan of the tube and In line with the preferred axis of vibra- 
tion. A flat circular piece of alimlmm was bonded to the member at mldspan 
for use as a target for the proximity probe. A photo llnagraph strip chart 
plotted the relative deflection versus time using deflection data froo the 
proximity probe. 


Evaluation of End Restraint Stiffness 

To evaluate the end restraint rotational stiffness K experimentally, the 
setup shown schematically In figure 4 was used. The base plate with the clevis- 
blade Joint (a) described earlier Is bolted to the test bed. A rigid round 
steel bar (b) Is screwed Into the top of the clevis. A cable (c) Is attached 
to the steel bar near Its top end and looped around the pull^s d, e, f, g, and 
connected to the connection blade close to the bottom, as shown In the figure. 
When a load W Is suspended through the pulley f. It Introduces equal and 
opposite forces In the horizontal portions (c and h) of the cable, thus 
cancelling the shear In the connection and the steel bar between the points of 
cable attachment while applying a pure bending mocsent to the clevis-blade 
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connection. Ihe ellealnatlon of the force Is necessary since the connection 
rotational spring stiffness K Is the slope of only the pure mooent versus the 
rotation relationship. A dial gage ms mounted at 12 In. from the base plate 
to measure the horizontal deflections of the steel bar b for successive Incre- 
ments In tlv^ applied load W. The deflections thus measured were used for 
computing the connection rotations. The maxImtiB value of U applied was 
2.6 lb. Five tests were conducted by attaching the shear- cancel ling cable to 
the lower clivis-blade attacheient screw, and another five were conducted by a 
similar attachment to the upper screw. The sane nuaber of tests were also 
conducted on each of the top and the bottom end fixtures. The average of these 
20 tests resulted In a linear moment- rotation relationship from which the 
following value of K was obtained: 

K » 1,974.3 1n-lb/rad1an (49) 

and was found to be similar to that obtained from approximate theoretical 
stiffness calculations using simple bending theory. 

Natural Vibration Tests 

A number of natural vibration tests were conducted on the partially 
restrained tubular member shown In figure 3 both with and without Internal 
passive damping devices. Figure 5 shows schematically the tubular member with 
(a) no Internal passive damping device; (b) with an Internal passive daupper In 
the form of a single lead shot weighing 0.47 ^ attached at midheight to a 
nylon fishing string weighing 0.17 gms which In turn Is anchored Internally to 
the fix ures at both ends of the member; (c) a mass- string system with a lead 
shot at each of the locations 0.25L, 0.50L, and 0.75L from the bottom of the 
member; and {A) a mass- string system with three lead shots at the member mid- 
span. Tho same nylon string was used In each arrangenent, and each Individual 
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lead shot In arrangetoents (c) and (d) weighed the same as the one In arrange- 
ment (b). The deflecti on-time plots were obtained by “plucking" the meaiber 
manually by a small amount either toward or away from the proximity probe at 
the midheight. Th' plot generated was used to calculate the frequency of the 
member, and the system damping. Although no Internal passive dan^jer was 
provided In arrangerent (a) the deflection-time relationship exhibited a 
gradual though slow decay due to the Inevitable presence of structural 
damping. This structural damping Is also present In arrangements (b), (c), and 
(d) besides that due to the Internal dampers. 

Figure 6 shows the cross-sectional view In "POSITION 1“ of the member at 
Its midspan relative to the origin 0 of the x,y,z coordinate system and the 
proximity probe. The deviation of the cross-sectional centroid from 0 Is due 
to the member Imperfection which diminishes to zero at the member ends. Since 
the direction of the applied Initial deflection (plucking) Is either along Oz 
or parallel to It, It was decided to study the effect on the member ctynaralc 
response of the various positions of the midspan cross section as shown In 
figure 6. The member had a near-con^plete restraint with respect to axial 
displacement at the ends owing to the rigid test bed at the bottom and a rather 
rigid bracket at its top. Thus, for example, If the member was plucked along 
Oz with position 2 as the initial location, some small amount of axial 
compression was Induced causing the member to vibrate outside of the plane Oz 
and also wobble to some extent. Although these three-dimensional effects are 
not accounted for in the analysis presented earlier In this memorandum, the 
response of the member parallel to or In line with Oz recorded experimentally 
cofl^ared favorably well with the corresponding analytical results as discussed 
later. The other Imperfection positions, 3 and 4, In which the member was 
tested are also shown In figure 6. 
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Table 1 presents a sumaary of the results based on convicting experlaents 
on the partially restrained meober with various loperfectlon positions (see 
fig. 6) and In the absence or presence of damping devices (see fig. 5). The 
test results are nunbered 1 through 8 each of which Is an average of three 
separate tests. Tests 1 through 4 were conducted In the absence of damping 
devices for the four 1n?>erfect1on positions 1 through 4, respectively. The 
aiQount of structural damping was measured by the logarithmic decrement method 
(see ref. 6) and Is sjnoarized 1n appendix D as It applies to the problem under 
consideration. The damping ratio c defined In appendix D Is also given In 
table 1 and Is calculated from the experimental deflection- tine response curve 
of the member at Its oldspan. The natural frequency of the member Is also 
calculated using the same curve and Is tabulated as fg. Tests 5 and 6 were 
conducted with the damping device (b) shown In figure 5 with Imperfection posi- 
tions 1 and 4, respectively. Tests 7 and 8 were conducted with Imperfection 
position 4, and with damping devices (c) and (d), respectively. 

Effect of Inqjerfectlon on Member Vibration 
In the absence of damping devices, the effect of the various Imperfection 
positions shown In figure 6 on the member dynamic response may be observed by 
coB^iarlng the c and the fg values for tests I through 4. The c values 
due to structural damping range from 0.0034 to 0.0183. The fg values vary 
from 3.16 to 4.10 Hz. It should also be noted that the smallest c value 
(0.0034) among these four tests does not correspond to the maxlraun fg value. 
It Is apparent, therefore, that the effect of the Imperfection position on the 
natural frequency of the raember Is significant. The axial end restraints, and 
three-dimensional deflection and wobbling of the member even though It Is 
plucked only In one plane, appear to be the main causes of the observed varia- 
tion In the fg values. The dynamic spatial response of the member brings 
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Into play the strong axis flexural stiffness of the end restraints as well as 
additional spatial damping effects which alter the fg values. In the 
presence of dafl^)1ng device (b), tests 5 and 6 show that the Imperfection posi- 
tion 4 results In higher damping and natural frequency values (e.g., larger 
C and fg values) as compared to those for the position 1. This t/P® 
behavior may again be attributable to the attendant spatial effects. 

Effectiveness of Passive DaJn^ers 

Test 4, 6, 7, and 8 were all conducted with imperfection position 4 (see 
table 1 and fig. 6) with damping types (a) through (d), respectively, as shown 
In figure 5. The ? and fg values are given In table 1 and the envelopes of 
the deflection-time response curves are presented In figure 7 In which a Is 
the normalized raldspan amplitude parallel to the Oz axis (see fig. 6) and t 
is the time In seconds. The outermost pair of curves In figure 7 represented 
by solid lines corresponds to damping type (a), that Is, structural damping 
only. The curves for daji^lng types (b), (d), and (c) are represented by 
dashed, dash-dot, and dash-dash-dot lines, respectively. Clearly, the damping 
type (c) with a lead shot at each of the locations x ° 0.25L, 0.50L, and 0.75L 
is the most effective of the ones considered in terras of Its ability to absorb 
energy and thereby reduce vibration. It should be noted, however, that the 
dSitping types (c) and (d) result in nearly Identical A-t envelopes up to 
about t = 1.5 sec, after which type (c) shows a slight superiority over type 
(d), although they may be considered practically the same. Furthermore, as 
seen from table 1, the c value for type (d) Is greater than that for type (c) 
while the fg value for type (c) Is smaller than that for type (d). Again, 
the spatial effects may be a cause of such response. 
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CQ»>ARISON OF THEORETICAL AMD EXPERIMEMTAL RESULTS 

With relatively small values of the damping ratio z such as those for 
tests 1 through 8 given In table 1, a significant variation of the natural 
damped frequency of the member with c cannot be expected as evident from 
equation (16). Thus, It Is concluded that the variation In fg values as seen 
from table 1 Is attributable primarily to the spatial effects mentioned 
earlier. To see how well the frequency predictions from the theoretical 
analyses presented earlier In this memorandin compare with those observed 
experimentally, the physical dimensions, the material properties, the end 
restraint stiffness, and the c values from table 1 were fed Into the computer 
programs developed specially for this purpose. The fg values from the exact 
solution, the approximate analytic solution (Involving eq. (28)), and the 
finite-difference solution for all the eight tests were found to be 3.58, 3.63, 
and 3.52 Hz, respectively. The theoretical analyses were conducted on a ajeober 
with no Initial Imperfections. These theoretical fg values are In vei 7 good 
agreement with the average experimental fg value of 3.42 Hz for the eight 
tests In table 1. Also, the experimental fg values are In the range 3.16 to 
4.10 Hz and all of the theoretical predictions are within this range. 

A comparison of normalized deflection- time curves for a time Interval of 
1 sec obtained from each of the three analyses was also made to those obtained 
experimentally. During the time Interval considered, the theoretical predic- 
tions were found to be In very good agreement with the experimental results. 

It should be recognized, however, that without Incorporating Into the 
theoretical analyses all of the existing spatial effects mentioned earlier, the 
theoretically predicted deflection- time curves will deviate from the 
experimental ones over longer Intervals of time (such as of the order of 2 to 3 
<^ec). That this would be the case may be seen by comparing the c values and 
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the A-t envelopes In figure 7 for tests 7 and 8. Any of the three 
theoretical analyses presented for a oes^er irlthout spatial effects trill 
predict a greater dec^y trfth time of the deflections for test 8 <^e to its 
larger c value of 0.0354 compared to that for test 7 trlth a c value of 

0.0296. The A-t curves for these two tests, in figure 7. show that bpyond 
t ° 1.5 sec approximately, the results are on the contrary. For a more 
accurate theoretical prediciton of the overall deflection- time response of the 
member, therefore, a three-dimensional analysis should be formulated. 

SUMMARY AND CONCLUSIONS 

The following is a brief susmary and some conclusions based on the outcome 
of the research reported herein: 

1. Explicit formulas for predicting the elastic buckling load and the 
natural frequency of colimns with partial rotational end restraints are 
developed. The buckling load formula gives results nearly identical to those 
from the exact, and the finite element techniques. The natural frequency 
formula also gives results nearly identical to those from the exact analysis. 

2. The natural frequency predictions from the formula developed, the exact 
analysis, and the finite- difference formulation presented are all in good 
agreement with those observed experimentally despite the presence of signifi- 
cant axial member end restraints and the attendant spatial effects in the 
member response. 

3. The member deflectiorhtime response, with a zero axial load, from the 
theoretical analyses is in ver 7 good agreement with that observed experimen- 
tally for about four cycles of vibration beyond which some deviation is 
expected due to the spatial and axial end restraint effects not included in the 
analyses. 

4. The three lead shot configurations provided considerably greater damping 
than the single lead shot damper. 
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APPEKOIX A 

nws IN THE MATRIX (ff EQUATIM (21) 

In this appendix, the 16 terms defining the matrix of equation (21) are 
given as follows: 


an - 

*13 * 1 


*12 “ 

«14 ■ » 


«21 “ 

- Elpf 


»22 “ 

- KSl 


«23 “ 

EIpJ 


«24 “ 

- 


«31 " 

COS p^L 


«32 “ 

sin PiL 


®33 “ 

cosh P 2 ^ 


«34 “ 

sinh 


a4i = 

EIpi cos p^L - 

Kp][ sin p^L 

«42 “ 

EIpi sin piL + 

Kp^ cos p^L 

®43 “ 

Elp| cosh p£L 

+ Kp 2 sinh pgL 

«44 “ 

EIp| sinh p^L 

Kp 2 cosh P 2 L 
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APPENDIX B 


TERMS OF MATRIX Cr] 

As an exas!p1e, the various terms of the niatrix [r3 used In e<^at1on (39) 
are defined In this appendix for H ° 6. The various teros are as follows: 



■ "55“ 

■ 6b^ - 

2b2- 

2bs - b4 - bi bg 

"12 

“ "23 

" "34 

" "45 

“ - 4bj + bg 

"l3 

""24 

""35 

" ^1 


"14 

" "15 

" "25 

B 0 


CM 

""33 

" "44 

“ 6bj^ 

- 2b2 - 2b3 . b^ 

"iJ 

""JI 






The terras b 3 and are defined by equations (42) and (43), respectively. 
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APFDCOIX C 


DERIVATION EXPRESSION 

The expression given by equation (45) Is derived In this appendix. 


With: 


k ["tN-*!) - •(‘(•“J] • ^ t'l-®’ *r^ (*r®) 

fit 


® bt^ ' 


suppose that equation (1) also holds on the Inititl 11ne» that Is, 


i-i) 


4 

EI^(x,.0) ♦ P 
dx 


Lex 


+ C^(x,.0). 0 


(C2) 


If ^[x^} In equation (6) exists, then: 

2 rr 4 bP 2 1 ^ 

ew/ Eibw, .. Pe«/ Cbw, 

2 " ■ p 4 ^*1’°^ p et 

8t 0x Lex j 


El d" f . P{/ , . 


(C3) 


dx 


dx‘ 


In which the Initial condition (7) Is already used. Using the central differ* 
ence expressions for the second and fourth order ordinary derivatives of 
^(x^) In equation (C3) results In: 

t*1'0) ‘ "VJ ^^1-2 ‘ * ^+2^ 

St n 


P 

p 


7 (* 1-1 ■ ^*1 * ♦ i » i > * *1 

.n 


(C4) 
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Substituting equation (C4) Into equation (Cl), and reerranglng gives: 


(♦(.J - ♦ «♦, - ♦ Vz) 


(At)^ P 
2 P 






(CS) 


which takes the fonn of equation (45) if condensed by using equations (46) 
through (48). 


APPEMOIX D 


EQUATION FOR DAWIKO RATIO 

The logarithBlc decreaent method (see ref. 6) tas used to calculate the 
damping ratio c defined as: 


C “ ^ (01) 

c 

In which C Is the damping coefficient and Is the critical damping factor 
given by: 


= 2 JeI^ (D2) 

The damping ratio was calculated using the following equation: 

C “ — (03) 
\/( 2 «)^ + 6 ^ 

In which 6 Is given by: 

*1 

6 ° in — (04) 

2 

where and X 2 are the successive aoplltudes of vibration as defined In 
reference 6. Each c value was obtained by taking an average of five 6 
values from each deflection- time plot. 
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